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A hypergeometric treatment to explain the nonlinear trne behavior of 
rednndant constraints on a straight elastic rod 


Giovanni Mingari Scarpello * Daniele Ritelli ^ 


Abstract 

In theory and practice of elastic straight rods, the statically indeterminate reactions acted by perfect 
constraints are commonly believed not to depend on the flexural stiffness EJ. We solve exactly two elastica 
problems in order to obtain hypergeometrically (helped by Lagrange, Lauricella, Appell), the true displace¬ 
ments upon which the forces method is founded. As a consequence, the above reactions are found to depend 
on stiffness: the presumptive independence credited as general, is far from being always true, but, quite the 
contrary, is valid only within a first-order approximation. 

Keyword: Non-Linear Rod Theory, Statical Redundancy, Elliptic and Hypergeometric Integrals, Lauricella 
Functions. 


1 Introduction 

In 1691, Jakob Bernoulli proposed to find out the deformed centerline (planar ’’elastica”) of a thin, homogeneous, 
straight and flexible rod under a force applied at its end. He established, Curvatura laminae elasticae (1694), 
the differential equation of the centerline, say y = y{x), of the rod ab appenso pondere curvata to be: 

, x^Ax 

His nephew, Daniel Bernoulli, computed the potential energy stored in a bent rod and, in 1742, wrote to L. 
Euler that the rod should attain such a shape as to minimize the functional of squared curvature. Accordingly, 
Euler dealt elastica as an isoperimetric problem of Calculus of Variations, De curvis elasticis (1744), arriving at 
elastica’s differential equation and identifying nine shapes of the curve, all equilibrium states of minimal energy. 

Anyway, C. Maclaurin had realized the elastica’s connection to elliptic integrals, see A treatise on fluxions 
(1742), § 927 entitled: ‘^The construction of the elastic curve, and of other figures, by the rectification of the 
conic sections”. In 1757, Euler authored a paper, Sur la force des colonnes, concerning the buckling of columns 
again, where the critical load is approached through a simplified expression to the elastica’s curvature. 

Some analytical solutions to elastic planar curves through elliptic integrals of the Hrst and second kind, 
can be read at m, a and [IH]. Almost half a century ago, [7], collected many problems on the base of the 
constraint type, and all solved by the elliptic integrals of first and second kind, while the third kind and Theta 
functions appear marginally for 3 — D deformations. In order to set the record straight, the concept of rod will 
be recalled from [5]: 

A “rod” is something of a hybrid, i. e. a mathematical curve made up of “material points”. It 
has no cross section, yet it has stiffness. It can weigh nothing or it can be heavy. We can twist 
it, bend it, stretch it and shake it, but we cannot break it, it is completely elastic. Of course, it 
is a mathematical object, and does not exist in the physical world; yet it finds wide application 
in structural (and biological) mechanics: columns, struts, cables, thread, and DNA have all been 
modeled with rod theory. 

In common practice, a simplified expression to the elastica’s curvature is currently used, while the real exact 
expression leads to nonlinearities. 
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The link between elliptic functions and elastica was always understood to be close, to induce the elliptic 
functions to be plotted through suitable curved rubber rods [5], while the paper [T] is founded on the use of 
the Weierstrass p function. Even in contemporary literaure the touchs of the special functions with elastica 
problems are very close, [?]■ In fact, the situations one can meet require either a numerical m, or an elliptic 
approach Elillilllll] or a hypergeometric m- In this paper we will follow the last one. 

2 Aim of the paper 

There is a large amount of literature providing the true elastica by solving the relevant nonlinear deflection 
equation (EH), where the rod’s curvature is forced by a bending moment. As a matter of fact in m , we 
solved the nonlinear deflection problem, as a pure strain analysis, for slender rods by means of Lauricella 
hypergeometric functions. 

Nevertheless, as far as we are concerned, the nonlinear, deflection analysis has never been used for a further 
static analysis until now. 

We refer to statically indeterminate structures, of which the degree of redundancy is the excess of constraints 
beyond the strictly necessary. This terminology is currently used as one can see in the literature, see for instance 
[U, or the Wikipedia entry “Statically indeterminate”. 

Such systems require additional equations and cannot be analyzed through the equilibrium equations alone. 
On the purpose, a frequent approach, the so-called forces method, obtains them by transforming the given 
redundant structure into a so-called simple one, i.e. a statically determinate one, by replacing the redun¬ 
dant constraints with the unknown reactions. This will not introduce alterations in stress distribution and in 
deflections. By decoupling such a multiple load condition over the simple structure in each load stage, each 
corresponding displacement is computed as a function of the system parameters (say E, J, L, ...) and the re¬ 
dundant unknowns. Finally, one shall require the consistency, namely that the displacement composition meets 
the real situation. And so we can have consistency to a zero displacement (perfect constraint), or to a certain 
imposed value (anelastic yield) or to a damper effect (elastic yield): in such a way the necessary elasticity 
equations are finally obtained. All the computing work one does about such deflections of straight rods of span 
L is grounded upon the small strains assumption: if y{x) is the deformed centerline’s equation, for each x it 
shall be: 

with 0 < X < L. (2.1) 

In such a way, all further steps -and the statically indeterminate reactions- lose their validity when the load 
intensity, beam slenderness, or high flexibility, causes the relevant strain-even kept in the linear elasticity (effect 
proportional to its cause) to be such that (EB becomes. In such cases, the exact centerline curvature is necessary 
and the nonlinear ODE m of deflections is often employed in advanced applications regarding the strength 
of materials in the contexts of aerospace or satellites, see [5]. 

We are just going, on the ground of these new considerations, to compute the statical indeterminate unknowns 
through a nonlinear deflection analysis. 

3 Statement of the problem 

Let us tackle a L-long, thin rod whose end is clamped and A free. We put at A the origin of a {x, y) cartesian 
reference frame with x along the undeformed rod, from the wall to A; and y downwards, normal at A to x, see 
Figured! Our basic assumptions are: 

Ai the rod is thin, initially straight, homogeneous, with a uniform cross section and uniform flexural stiffness 
EJ, where E is the Young modulus, and J, the cross section moment of inertia about a “neutral” axis 
normal to the plane of bending and passing through the central line; 

A 2 the slender rod is always charged by coplanar dead (not follower) loads; 

A 3 Linear constitutive load holds: so that the induced curvature is proportional via 1/{EJ) to the bending 
moment. 

A 4 The shear transverse deformation is ignored. 

A 5 A stationary strain held, by isostatic equilibrium of active loads and reactive forces takes place, and, due 
to rest of static equilibrium, no rod element undergoes acceleration. 
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The key Bernoulli-Euler linear constitutive law connects flexural stiffness EJ and bending moment ^ to 
the consequent curvature x of the deflected centerline: ±EJx = ^{x). With our (traditional) reference frame 
{x,y), the curvature always bears the opposite sign to the bending moment: 


EJ -= -^{x) (3.1) 

where y = y{x) is the unknown elastica’s equation and 

_ ^ ^ 

Ux — 1 ■) Vxx — 1 O • 

ax ax^ 

To ()3.1|) a useful dress can be given, so that it shall be minded that the Cauchy problem: 


A [ I 

\ \/l + vl) 

y{L) = 0 


yx{L) = 0 


— ^—^{x), 0<x<L 

EJ 


(3.2) 


can be solved in such a way. Putting: 


the solution of (EJ) is: 


y{x) = - 


H{x) := f .^(C)de 

J X 

H{0 


de 


jEjf-H^O 

In order to have a solution defined for x € [0, T] we have to require that 

- 1 < Eh(:,) < 1 

which is a general prescription on the load effects and will be defined in single cases. 


(3.3) 


(3.4) 


3.1 Nonlinear deflection, a first example 

For a, L, E, J cantilever tip-sheared by P, the approximate linear theory provides a tip’s vertical displacement 
i5o = {PL^)/{3EJ). In such a case, in (13.21) we put ^{x) = —Px and defining the nondimensional quantities: 


PL^ y ^ X 

^ ^ 3^7’ ^ ^ L’ ^ ^ L 


it becomes Sq = 2y/3. 

Formula (2.32) of our previous paper [TS] for the exact free-end displacement provides: 




3(1 - 


e 



where d(/i) is a complicated function involving elliptic integrals of first and third kind according to formula 
(2.29) of the same |15j . They can be expanded with respect to fi, with initial point /xq = 0; then, recalling the 
Maclaurin first-order expansion with respect to /x 




3(1-a^2)3/2 




0{y^) 


(3.5) 


one obtains eventually 

m = ^{2-3^ + e)+0 (3.6) 

i.e. the usual approximation which for ^ = 0 (free tip) provides the same deflection (5o as above. 
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Figure 1: Cantilever’s kinds of deformed centerline and transverse displacement A 



Figure 2: The heavy cantilever supported by a roller (a), the relevant simple structure (b) and the loads 
decoupling (c) and (d). 


Figure [T] shows qualitatively the behavioral differences, for a L-rod OA under a dead load P at the tip, 
between its typical approximate elastica OA and the true elastica OA'. As a matter of fact, the exact treatment 
highlights a horizontal displacement of the free tip which is given by 

A = L — inf [ 

0<y<a;(L) 

and which has been computed (sect. 2.1 and 3.) in [T5] . 

What we are now going to do is to employ a displacement computation like this for gaining exact elasticity 
equations to solve statically indeterminate systems. 

4 The statically indeterminate heavy cantilever supported by a roller 

Our first statically indeterminate system is a rod (Figure [5] (a)) under a uniform load q, whose one side A is 
built-in, while the free tip B is supported by a roller: for such a 1-redundant rod we assume as unknown the 
reaction ^ just acted by the roller. Then the statically “simple” rod is of (Figure [2] (b)), whose analysis can 
be decoupled according to Figure [5] (c) and (d) and tackled separately. 
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4.1 First sub-system: the heavy cantilever 

Let us consider first the subsystem of Figure [5] (c). For the bending moment we have ^{x) = —{q/2)x^. 
Specializing p.4l) for this load, it will imply q < (SEJ/L^. Thus (13.31) read^ as 


y{x) = f 


q{L^ - y3) 


: dy = 


2EJ 




qL^ 




2/3 


du 


TiL^-y^y 

Now using the Lauricella functions’ integral representation theorem, we infer the integration formula 

a, —a, ab 


(4.1) 


Jo (1-u2)1/2(1_6i()2/3^^ 2^^ 

Clearly (14.21) is the key to evaluate the integral in (14.11) . so that 


a-^( 3 ) / 2; 1/2,1/2, 2/3 


(4.2) 




2; 1/2,1/2, 2/3 




dHb) 


which provides the deflection law for each x. We are interested to take a: = 0 in (I4.1b|) getting 

L^q L^q \ 


36EJ ^ I a 


6EJ’ 6EJ' 


(Hil l 


Notice that for a; = 0 the third argument of Lauricella’s becomes 1, so that such a function collapses into an 
Appell one and from there it can be reduced to a 3 F 2 hypergeometric function, as shown in section | 6 l In such 
a way we will be allowed to describe the strain by means of a function of only one variable. We can simplify 
(I4.1bl) using the reduction formulae (16.31) and (16.41) in sequence, getting 




!,1 


2 ’ 2 

1 5 
6 ’ 3 




m{Ejy 


dHli) 


which provides as a function of q, L, E, J the free tip true deflection of rod of Figure [2] (c). 

4.2 Second sub-system: the tip-sheared cantilever 

Allow us consider the subsystem of Figure [5] (d). The bending moment being ^(x) = S^x, so that (13.3p reads 
as 


y{x) = / 

J X 


) , 3; r _ u _ 

^(2EJY-S-Hy^-L-Y * + " 


With the integral representation theorem we also infer the integration formula 

' 2 ; 1 / 2 , 1 / 2 , 1/2 


a^„(3) 


/ _ci?/ 

/o (1-62 u2)1/2(1+cu)1/2 2 ^ 


ab, —ab, —ac 


(4.3) 


(4.4) 


Clearly, (14.41) is the key to evaluate the integral in (14.3p . giving 


y{x) = 


8EJL 


^ D 


2 ; 1 / 2 , 1 / 2 , 1/2 

3 


{x^-L^)^ 

2EJ 


(x^-L^)^ 

2EJ 



6^3) 


^The approach of computing the exact tip position of the rod would be extremely long and capable of hiding the sense of this 
research deviating by its central purpose. In such a way the (very small) horizontal deflection of the tip has been put to zero. 
Anyway, an example where the rod is assumed inextensible and the axial deformation is kept in account, can be seen in our previous 
paper [14] sect. 3, formula (3.12). 
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We are interested to take x = 0 in 1130, getting 


y(o) = 


8EJ ^ 


2 ; 1 / 2 , 1 / 2 , 1/2 

3 


2EJ ’ 2EJ 


We can simplify (I4.3hll using the reduction formulae (16.3|) and (16.41) in sequence, getting 


y(o) 


' 3EJ 


3-r 2 


I, I <Ejy) 


(03;) 


6 JH) 


4.3 Consistency 

In order to evaluate the redundant unknown notice that the addition of two displacements given by 030 
and (03E1) shall be zero, arriving at the transcendental equation: 


3Lq 3 F 2 


1 1 3 

2 ’ ’ 2 

I 5 

6’ 3 




36{Ejy 


-SJTgFa 


1 1 3 

2 ’ ’ 2 
5 I 

4, 4 


4{Ejy 


= 0 


(4.5) 


in our unknown reaction ^. Of course, (H3|) cannot be solved analytically, but numerically. Nevertheless we 
will provide an analytical approximate solution. In fact, if we approximate both hypergeometrical series only 
by their first term which is equal to 1, the solution to (14. 5|) is approximately given by: 



according to the linearized theory. 

For a higher approximation, we can proceed through the Lagrange series reversion^. 
For the purpose, in (14.51) we put: 


jr = 


2EJ 




so that (H31) can be written as 


/(iT) :=(r3F2 


/ 1 1 3 

( 2 ’F 2 


ill 

2 ’ ’ 2 

1 ^ ^ 

\ 4’ 4 

j 16 EJ 

7 5 

6’ 3 


6„2 


L'^q 


MEJf 


= 


0313) 


Function /:] — !, 1[—>■ R appearing at left hand side of (I4.5b l is odd, bijective and then invertible. Notice that, 
due to the hypergeometric series divergent behaviour at the edge of the convergence disk, we have: > +00 

as —>■ 1 “. 

In Figure|3]we depict continuously the plot of /(i^) and as a dotted line the inverse By Mathematica, 

based on the Lagrange inversion theorem, we can obtain the series expansion with respect to 2^ of the root of: 
/(i^) = . We give hereinafter the relevant computer instructions in Figure S] 

After this, we will expand the right hand side of (I4.5b l. evaluate the approximate solution and expand it in 
power series: 

Now, going back to our unknown ^, we can get an approximate solution of (I4.5b l: 


8Lq 153L^g3 47277L^^q^ 4596133617 L^^q'^ 

~ ~ 143360(^7)2 “ 267177164800(£;J)4 “ 200130658356428800(F;J)6 ’ 


The plot in Figure [5] shows the behavior of ^ versus the number of terms of (|4.6D wherefrom the linearized 
solution is coming for n = 0 and overestimates the solution; the addition of further terms provides us with a 
better approximation which is stable just after n = 7. 

^Lagrange, J.L. (1770): ’’Nouvelle methode pour resoudre les equations litterales par le moyen des series,” Memoires de 
TAcademie Royale des Sciences et Belles-Lettres de Berlin, vol. 24, pages 251-326, a paper of 1766 but issued four years later. The 
Lagrange inversion theorem used here gives the Taylor series expansion of the inverse function of an analytic function; it shall not 
be confused with Lagrange reversion theorem which provides formal power series expansions of certain functions defined implicitly. 
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Figure 3: /(^) and 


ini3]:= s = Series l^x HypergeometricPFQ 1^11, —0/ ^0} 
tt = Normal[inverseSeries[s, y] ] 

27y3 99 873 y5 781 224 831 y'^ 620 107 856 469 y^ 

Outi;4]= y - -+--+ - 

70 509 600 7 455 448 000 10 855 132 288 000 


Figure 4: Series for 


2 Q 3 q 3 , 7 5 , L‘ 

in[5]:= solappr = - tt / . y-> - HypergeometricPFQN —, 1, — - 

l2 16 Q >-^2 2-J ^6 3-^ 36q2 

Normal[Series[solappr, {q, 0, 7}]] 

3Lq 4 596 133 617 q’ 47 277 l“ q= 153 l’ q^ 

Oul(6]=---- 

8 200 130 658 356 428 800 Q® 267 177 164 800 143 360 0^ 


Figure 5: Approximate solution of (|4.5b ) 


375 - • 

370 - 

365 - 

360 - 

355 - 

350 - 

. 

5 10 15 20 25 

Figure 6: The statically indeterminate heavy cantilever supported by a roller with L = 1, q = 1000, EJ = 200, 
^ values versus number of terms of its expansion 
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Figure 7: The doubly built-in cantilever (a), the relevant simple structure (b) and the loads decoupling (c) and 
(d). 


5 The statically indeterminate two-side built-in heavy rod 

Our second problem, Figure [71 will concern the constraints’ redundancy of a doubly built-in slender rod AB 
subjected to a uniform continuous load q. Such a system, for a general load, is three times redundant; for 
a purely vertical but non-symmetrical one, two times. In our case the actual degree will be only one. We 
assume as unknown ^ the bending moment acted by one of the built-in supports (see FigurejT] (b)). The same 
constraint will provide a vertical reaction equilibrating just half a load and then of value (?T/2. In such a way, 
the bending moment at the cc-section will be given by: 

... , qLx ^ 

^{x) = —^- 0 < X < L. 

As before, we are going to analyze apart the subsystems (c) and (d). 


5.1 First sub-system: a heavy tip-sheared cantilever 


Let 


^{x) = 


Lqx 


qx^ 


for the rod of FigurejT] (c). Observe that (13.41) reads as 

-1 < ~ + 2a;) < 1 

which will imply the condition q < VIEJjL? to be met. 

Accordingly, (13.31) lead^ to the hyperelliptic integral: 



_ <l{L-y)^{L + 2y) _ 

v'144(F;j)2-g2(L-j/)4(L + 2y)2 ^ 


(5.1) 


and in particular, for the free tip, obviously: 

q{L-y)^{L + 2y) 


2 /( 0 ) = - / 

Jo 


0 v'144(AJ)2 - q^{L - 2/)4(L + 2y)^ 


dy := -I{L,q,E,J) 


(I5.1b l 


^Being the double clamping made for guarantee two fixed positions, the maximum range of integration is from 0 to L. Then we 
fix X € [0, L], which implies the deformed rod is extensible under its load. 


8 



















































100 



Figure 8: The doubly built-in cantilever with L = \, q = 1000, EJ = 200: ^ values versus number of terms of 
its expansion 


5.2 Second sub-system: a cantilever under a bending moment at its free end 

Let ^{x) = ^ for the rod of Figure[7] (d). In such a way (13.31) is solved by elementary methods: 


, , JiEjy - Lf - EJ 

y{x) = -- 

which gives in particular 

y/{EJY - _ EJ 

yiO) = - ^ - 

and then we succeed in evaluating the displacement provided by the redundant unknown ^ if it would act 
alone. 



5.3 Consistency 


To compute the redundant unknown our constraint being perfect, we shall require again the addition of 
two displacements given by (IS.lbp and (ICT) shall be zero. The computation of the integral dSIH), because of 
the cubic, would lead to Lauricella higher order functions. By expanding in powers the integrand in (IS.lbl) , we 
get the approximation: 


I{L,q,E,J) ~ 


24:EJ 


1 2 
2 ’ 3 
5 
3 


\ 

m{EJY j 


(5.3) 


Then, solving to the displacement equation: 


^/{EJf - - EJ 


= I{L,q,E,J) 


2EJT{L,q,E,J) 

L^+I^{L,q,E,J) 


Putting to I{L, q, E, J) its approximation given by (15.3p and after a further development, we get: 


L^q ^ ^ TTL^'^q 

U ^ 34560(F;J)2 19906560(^7)4 ‘ ‘ ’ 


which holds even powers of EJ. Figure [5] shows the behavior of such a solution with the number of terms of the 
series (EUl): the linearized solution is that for n = 0 and underestimates the 3E ; the addition of further terms 
provides us with a better approximation which is stable after n = 11 . 


6 Conclusions 

In order to get a better clarity, we will refer to a sample problem. Let us take into account a slender rod of 
length L = Im and a square uniform section of side s = I0“^m, so that J = s4/I2 = The material is 

steel with E = 2.1 ■ lO^^N/rn^, so that each section has a flexural stiffness EJ ~ 200Nm^. 
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Let us consider first the statically indeterminate unknown relevant to Figure [2]: according to the approximate 
linear theory, it is given by 3gL/8, namely 375iV for our sample problem. We take a unit continuous load of 
q = IQ^N/m which is compliant to the prescription of being less than 6EJ/L^. 

Our nonlinear approximate solution is provided by a relationship of kind (14.61) , namely: 

°° r6n+l„2n+l 

^ - Va - _ ^ _ 

n=0 ' ^ 

The reaction ^ of the statically indeterminate roller is then computed via a power series involving length, flex¬ 
ural rigidity, unit load and the relevant coefficients, some of which are shown by (IMl). Therefore Figure [5] shows 
that the linearized approach overestimates about 7.8% the “reference value” provided by the (less approximate) 
nonlinear treatment which converges from n = 7 on. 

Passing to the second statically indeterminate system (see Figure [7]) we took the same unit continuous 
load q = IQ^N/m which is compliant to the new prescription of being less than \2EJ/L^. According to the 
approximate linear theory, the statically indeterminate bending moment (hereinafter called again ^ ) exerted 
by one of the building-in constraints, is traditionally accepted as qL'^jl2, namely 83.33iVm for our sample 
problem. But a proper number of terms of (15.4p 


r6n+2„2n+l 

r = ^ 

2^ " (F;j)2n ’ 

provides a value of 95.16 A^to which is 14.19% greater. 

The general conclusions suggested by all of the above are the following: 

Cl In both problems analyzed in this paper, the traditional linearized approach leads to significant errors on 
the statically unknowns. Such deviations increase for rods that are progressively more and more slender 
and flexible, see Figure HI Figure HI 

C 2 The above mentioned effects obviously lead to a remarkable impact on the axial bending stresses, say cr 
s. E.g. for the rod in Figured the highest bending moment is given by 37“^I2q and in that section the 
relevant edge tr is 16% lower than that from the linearized approach. 

C 3 For instance, referring to the rod in Figure [71 in the section of the highest moment, the relevant edge a is 
6 % greater than that from the linearized approach. 

C 4 It is commonly believed and written that the redundant unknowns do not depend on the flexural stiffness, 
but this is, generally speaking, not true. 

In fact, our final formulae dm and dm show undoubtlessly that this is due to having taken n = 0 , namely 
only the first term of the expansion whose subsequent terms hold the reciprocal of powers of EJ. Therefore the 
Belluzzi sentence, see [3] p. 363: 

Si osservi che quando i vincoli sono rigidi la reazione A non dipende da EJ, mentre ne dipende 
quando sono cedevoliQ 

can be accepted solely within a first order approximation. 

The linearized approach does its approximation at the very beginning of the treatment, so that leads to a less 
refined and elementary conclusion. On the contrary, we keep the nonlinear nature of the problem, generating 
a true solution which -even if we decide/need to approximate it at final stages- reveals new mathematical links 
and conclusions. 


Appendix: Hypergeometric identities 

We recall hereinafter something about the non too known Lauricella hypergeometric functions. 

“^Notice that when the constraints are rigid, the (redundant) reaction does not depend upon EJ^ whilst it does when they 
undergo some yield. 
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The first hypergeometric series appeared in the Wallis’s Arithmetica infinitorum (1656): 


a-b a • (a + 1) • 6 • (6 + 1) 

,F, (a. I.; c; *) = 1 + — ,T + ^ 

for |a;| < 1 and real parameters a, b, c. The product of n factors: 

(A)n = A (A + 1) • • • (A + n — 1), 

called Pochhammer symbol (or truncated factorial ) allows to write 2 F 1 as: 

(a)„ ib)n 


2 Fi{a,b;c-,x) = ^ 


(c) n! 

n=o '■ 


A meaningful contribution on various 2 F 1 topics is ascribed to Eulei@; but he does not seem to have known 
the integral representation: 


2 Fi{a,b;c;x) = 


r(c) 


^(1 — u) 


c—a—1 


r(a)r(c—a) Jq {1 — xu)^ 


■ du, 


really due to A. M. Legendr^l. The above integral relationship is true if c > a > 0 and for |a;| < 1, even if this 
limitation can be discarded thanks to the analytic continuation. 

Furthermore, many functions have been introduced in 19*’*' century for generalizing the hypergeometric 
functions to multiple variables. We recall the Appell Fi two-variable hypergeometric series, defined as: 


Fi 


a; 61,62 


00 00 


\ (a)mi+m2 (^l)iTii (^2)m2 ^ ^2 \ ^ T \ 1^1 

xi,x2 = 2^ 2^ - — -l*^2| < i 


■mi=0 7712—0 


(c) 


mi+m 2 


mil m 2 ! 


The analytic continuation of Appell’s function on C \ [1, oo) x C \ [1, oo) comes from its integral representation 
theorem: if Rea > 0, Re(c — a) > 0: 


a; 6 i, 62 


Xl,X2 = 


r(c) (1- 


r(a)r(c — a) Jq (1 _ xi (1 — X 2 u) 


f )2 


du. 


( 6 . 1 ) 


For us, the functions introduced and investigated by G. Lauricella (1893) and S. Saran (1954), are of prevailing 
interest; and among them the hypergeometric function of u G N+ variables (and n + 2 parameters), see 
m and m. defined as: 


F^^ (a, 6 i,..., 6 „;c;xi,...,a:„) := ^ 


- \-7nn{^l)7ni ' ‘ ’ (^n) 


(c)miH- \-7nn^l^- ' ' ' ^n! 


'^ri mi ^ ^ ^ mrr 

ct/l 


with the hypergeometric series usual convergence requirements |xi| < l,...,|xn| < 1 - 
If Rec > Rea > 0 , the relevant Integral Representation Theorem (IRT) provides: 


(a, 6 i,..., 6 „;c;xi,...,x„) = 


r(c) 


*(i — uY 


r(o) r(c — a) Jq (1 — xiu)*’i • • • (1 ~ XnuY' 


■ du 


allowing the analytic continuation to C” deprived of the cartesian n-dimensional product of the interval ] 1 , oo[ 
with itself. 

Finally, we provide here some reduction formulae employed throughout this paper. 


a; b, b, 
u T 1 


X, -X 1 = 2 F 1 


\;b, 


- +1 
2 ^ ^ 


( 6 . 2 ) 


^We quote three works: a) De progressionibus transcendentibus, Op. omnia, S.l, vol. 28; b) De curva hypergeometrica Op. 
omnia, S.l, vol. 16; c) Institutiones Calculi integralis, 1769, vol. II 

®A. M. Legendre, Exercices de calcul integral, II, quatrieme part, sect. 2, Paris 1811 
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Identity (16.211 follows from the integral representation theorem. In fact, if we return from the symbolic expression 
for the Appell function to its meaning in terms of integral representation, we see that 


a; b, b 
n -j- 1 


r(a + l) 

r(a)r(i) X {i-x^u^y/^ 


du 


r ^ j F I 

= ®/ T, - 9 M /9 = 2 F 1 

Jo (1 - x^vy/^ \ t +1 


Next, we consider some generalizations of the well-known Gauss summation formula: if c > a + 6 then 


2 F 1 


1 : b 


1 = 


r(c)r(c — a — b) 
r(c — a)r(c — b) 


moreover, we also have 


Fi 


a; 61,62 


l,x \ = 2 F 1 


a; 61 


1 2 F 1 


o; 62 
c — bi 


which, when c > a + bi gives, due to the Gauss summation theorem. 


a; 61,62 


\ _ r(c)r(c-a-6i) / a-M 

’ r(c-a)r(c-6i) " M c -61 


In similar way, it is possible to recognize the following formula, provided that c > a + 63 , which we used to 
simplify equation (OB . 


.(3)/ a;&F&2,63 


D 


]_r(c)r(c-a-^ ( a; 6 i ,62 

j r(c-a)r(c-63) M c-63 


x,y 


We also recall the reduction formulcQ: 


Fi 


a; 6 , 6 , 


x,-x\ = 3F2 


Q+1 a t 
2 ’ 2 ’ ^ 

c+1 c 
2 ’ 2 


(6.3) 


(6.4) 
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